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Abstract 


We present a NLO perturbative analysis of all available data on the polarized structure function 
gi(x,Q 2 ) with the aim of making a quantitative test of the validity of the Bjorken sum rule, 
of measuring a s , and of deriving helicity fractions. We take particular care over the small x 
extrapolation, since it is now known that Regge behaviour is unreliable at perturbative scales. 
For fixed a s we find that if all the most recent data are included Qa = 1.19±0.09, confirming the 
Bjorken sum rule at the 8% level. We further show that the value of a s is now reasonably well 
constrained by scaling violations in the structure function data, despite the fact that it cannot 
yet be reliably fixed by the value of the Bjorken sum: our final result is a s (mz ) = 0.120 j/yoog- 
We also confirm earlier indications of a sizeable positive gluon polarization in the nucleon. 
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Much experimental and theoretical work has been devoted in recent years to polarized 
deep inelastic scattering [|1|. Reasonably precise data on the polarized structure functions of 
proton i-i and deuteron ||-[@ have been collected down to values of x near and below 
x = 0.01 for Q 1 2 > 1 GeV 2 . Very recently precise results on the neutron structure function 
from scattering on 3 He targets have also become available || [Hj. The calculation of the 


kernels for the perturbative QCD evolution of polarized parton distributions has recently been 
completed to next-to-leading order (NLO) [[Dj, thus reaching the same level of accuracy as in 
the unpolarized case. Experience obtained from the small x behaviour of unpolarized structure 
functions observed at HERA JT2| is now sufficient to indicate at least qualitatively the sort 
of behaviour we might expect for polarized structure functions at small x. The purpose of 
this paper is to take advantage of all this accumulated knowledge and experience in order to 
extract from the data on polarized structure functions the polarized parton densities and their 
moments for comparison with theoretical expectations. In particular we will show that when 
all the data are included it is now possible to make a reliable test of the Bjorken sum rule 
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We will also investigate to what extent one can use the data to make an accurate determination 
of a s . 

For an experimental verification of the Bjorken sum rule one has to extract from the data 
the first moment of the difference of polarized up and down quark densities at some convenient 
value of Q 2 . Data taken at all kinematically accessible values of x and Q 2 . and on all available 
targets, contain information relevant for the reconstruction of polarized parton densities at a 
given Q 2 and ought therefore to be included. The complete NLO evolution kernels |JTT[ can be 
used to reduce to the same Q 2 data measured at different Q 2 for each £.[] Since the evolution 
equations [|T3j for partons at a given x and Q 2 depend only on the values of the parton densities 
at larger values of x and the same Q 2 , the necessary correction can only be performed through 
a general fit to all the data, which yields a set of polarized parton densities obeying the correct 
evolution equations |16], O] . However in order to perform a fit one must start with a particular 


ansatz for the parton densities at some reference Q q. Clearly the results of the fit will depend 
to some extent on the starting ansatz one adopts, and this dependence will induce an error in 
the computed first moments, and in particular in the Bjorken sum. Here we will devote special 
attention to this issue. 


Once the data are reduced to a common Q 2 for all x values, an extrapolation to unmeasured 
values at small and large x is needed in order to obtain the first moment. The extrapolation 
at small x is especially important |TB|.0 In most of the existing analyses, including those in 


the experimental papers, it has been performed by assuming a simple power behaviour based 
on Regge theory JllJ. This leads to a rather small contribution to first moments from the 


small x region, since the expected extrapolation is at most flat. But two main considerations 
are now severely undermining these attempts. First of all the data at small x for cq indicate 


1 In the past this has been done by assuming that the polarization asymmetries are Q 2 independent [m|, 
but this approximation is questionable given the current precision of the data. 

2 Note that the behaviour at small x of the input ansatz for the parton densities at Q fj is not relevant for the 

evolution correction, which only depends on x values larger than the smallest measured one. On the contrary 
the integration at small x that completes a given moment is very much dependent on the small x behaviour of 
the input distributions, as we shall see. 
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a clear departure from the nearly constant behaviour expected from Regge theory, pointing 
to a more singular behaviour which is moreover different for protons and neutrons. Second, 
we have learnt from HERA experiments that the rise of unpolarized structure functions is 
also much steeper than that predicted by Regge theory 0 - Although Regge theory seems 
to work reasonably well at scales typical of soft hadronic physics, at larger scales the effects 
of perturbative QCD evolution are superimposed to give a rise which gets steeper and steeper 


as Q 2 


increases 


A similar phenomenon is to be expected for the polarized structure 


functions ||l6l [22j| , so the results of a naive Regge extrapolation are not justified if one wants to 
consider Erst moments in the perturbative region: small x contributions to first moments can 
be relatively large, especially as Q 2 


'* 2 increases. 


Here we will discuss alternative extrapolation procedures and the errors associated with 
them. Our guiding principles will be the validity of Regge predictions at low Q 2 , the buildup 
with Q 2 of the effect of the QCD evolution and the limits on the size of polarized densities 
imposed by the corresponding unpolarized densities. From these starting points we will esti¬ 
mate the uncertainty in the small x extrapolation, which when combined with the evolution 
corrections and the more standard sources of error will allow us to quantify the extent to which 
the Bjorken sum rule may be tested using existing data. In practice we will do this by deriving 
from the data the value of (]a and the associated error for an appropriate range of values of a s . 

We will then consider the determination of a s from the polarized deep inelastic scattering 
data. Previous attempts in this direction [^3|] have assumed the validity of the Bjorken sum 
rule, and used a value for the Bjorken integral obtained from the first moments given by the 
various experimental collaborations, and thus based on naive Regge extrapolation at small x. 
However when the effects of perturbative evolution on the small x extrapolation are properly 
taken into account, the evaluations of the first moments must be revised: their errors then turn 
out to be considerably increased, and the determination of a s from the Bjorken sum rule no 
longer works so well. However, we are able to show that a much better determination of a s may 
be obtained if all the available data and not only the Bjorken integral are used in the analysis: 
the comparison of the data at small and large Q 2 in the measured range of x then leads to a 
reasonably precise measurement of a s . 

A similar procedure to that described for the Bjorken sum rule clearly also applies to the 
determination of other Erst moments, and in particular to the hclicity fractions carried by the 
gluon density and diEerent flavours of precisely defined quark densities. Such determinations 
using NLO evolution kernels have already been presented in refs. 

mm 


17, M. What is new here, 


beyond the inclusion of the recent new experimental results 


is a careful analysis of the 


dependence of the results on the assumed form of the input parton densities and thereby on 
the extrapolation to small x. As a result we confirm the conclusion, first reported in ref. |17| , 


that the data indicate the presence of a sizeable positive gluon polarization in the proton, as 
conjectured several years previously |25|-|27| . 


Finally, we explicitly discuss the possibility of a direct determination of a s through the 
Bjorken sum rule from the measurement of the isotriplet first moment of g\. If it were possible 
to determine this first moment very accurately, one might then be able to take advantage |23[ 
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of the fact that radiative corrections to the nonsinglet first moment are know to NNNLO, that 
is to order a^. We show that, with current data, even though the central value of a s obtained 
with this procedure is quite compatible with that of the NLO analysis the experimental error 
is considerably larger, due to the uncertainty in the contribution to the isotriplet first moment 
from the unmeasured small x region. 


1 Polarized Structure Functions and Partons 


We begin by summarising various results on the relation between structure functions and po¬ 
larized parton distributions, and their behaviour at small x, which will be important for the 
following discussion. 


1.1 Defining Polarized Parton Densities 


The structure function g\ is related to the polarized quark and gluon distributions by |17 


gfix, Q 2 ) = C NS ® A q NS + C s 0 AS + 2 n f C g 0 A g], (1) 

where (e 2 ) = nj 1 e 2 , 0 denotes convolution with respect to x, and the nonsinglet and 
singlet quark distributions are defined as 


A qNS = Xj 
2=1 


■py 


2=1 


( 2 ) 


where Ag* and Aqi are the quark and antiquark distributions of flavor i and Ag is the polarized 
gluon distribution. The evolution equations for the polarized parton densities are given by 


dt 


A qNS 


d_ ( AS \ 
Jt\Ag) 


^A-p, n , s ® A ?™> 

“•« (Pi 2 n,P°\ 

pi) 



(3) 


where t = logQ 2 /A 2 . The coefficient functions C and the polarized splitting functions P are 
now known at LO and NLO ||TT[ . Moments of coefficient functions and parton densities 
are defined as f(N ) = /q 1 dxx N ~ l f(x) and denoted by C(N, a s ), Aq NS (N , Q 2 ), ATj(N, Q 2 ) and 
Ag{N,Q 2 ). 

As is well known |l]], the dehnition of the singlet quark density AS(i, Q 2 ) must be carefully 
specihed, since its first moment is maximally scheme dependent. This follows from the fact 
that, due to the axial anomaly, a s (t)Ag(l, Q 2 ) is scale independent at LO. This implies that 
even asymptotically the ambiguity in A£(l, Q 2 ) is necessarily of the same order as its size. For 
a sensible comparison with the constituent quark spin one must thus define AS(1,Q 2 ) in such 
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a way that it is scale independent [25, 28|: AE(l,g 2 ) = AE(1). This is the definition we will 
adopt here. We then have 


IMQ 2 ) = [ dxgi(x,Q 2 ) = ^-[C NS (l,a s (t))Aq NS (l) + C s (l,a s (t))a 0 (Q 2 )], (4) 

Jo 

with a 0 the singlet axial charge: 

a 0 (Q 2 ) = AE(1, Q 2 ) - n f *®Ag(l, Q 2 ). (5) 

The higher moments of the singlet quark distribution are also scheme dependent, although in 
a less dramatic way. Various schemes were discussed in ref. JH], ^9] and the dependence of the 
results of the analysis on the choice of scheme was studied. Here we do not come back to this 
issue but instead simply adopt the AB scheme as defined in ref. [171 . 


1.2 Small x Behaviour 


In view of the need to extrapolate the data to x = 0 in order to compute moments, it is impor¬ 
tant to summarise the current understanding of the small x behaviour of structure functions. 
For the unpolarized singlet quark and gluon distributions the QCD evolution equations (|3|) lead 
to the following asymptotic behaviour at small x [| 


21 


xg 

xE 


-l/2 e 2 7CT -5C( 1 + £" =1 e y +1 c4), 


cr 


1 + E?=1 qp^cc 


r ^+1 / 


( 6 ) 


where £ = logx 0 /x, C — log (a s (Ql)/a s (Q 2 )), a = \/££, p = and the e terms indicate 

corrections from the n-th perturbative order (with n — 1 corresponding to NLO). It follows 
that the structure functions xFi and F 2 rise at small x more and more steeply as Q 2 increases, 
though, for all finite n, never as steeply as a power of x. For all other parton distributions / 
(/ — Qns, Aq N s, AE, A g) one has similarly [|2], |l(| 


/ - a - 1/2 e 2 7/-5/C(i + E" =1 e}p 2i+1 y). 


(7) 


Thus these distributions are less singular by a factor of x than the singlet unpolarized distri¬ 
butions eq. ([|), while the higher order corrections are more important at small x since the 
exponent * + 1 is replaced by 2i + 1; this is because the leading small N contributions to 
the anomalous dimensions at order cc* +1 are («s/( N — l)) 4 in the unpolarized singlet case, but 
N (a s /N 2 y for the nonsinglet and polarized distributions. 

The limiting behaviour (jo].[7;) implied by the evolution equations (J5|) at finite order in per¬ 
turbation theory would be spoiled if the series of higher order powers of logl/x were summed 
to all orders to give a powerlike behaviour in x, which would then overwhelm the leading terms. 
As is well known |F|, in the unpolarized singlet channel one may obtain a result as singular 
as x~ x , with A ~ 1/2, for xE and xg by summation of higher order singularities in the Regge 
limit of x —> 0 at fixed Q 2 (hence fixed a s ). However the meaning and the value of a fixed 
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a s are quite ambiguous, and it is not at all necessary a priori that such a singular behaviour 
is of relevance in the measured HERA region. In fact the experimental results from HERA 
show no evidence at all for this behaviour [2T|, [12]]. In principle the higher order terms could 
be more important for the nonsinglet and polarized distributions due to the 2i + 1 exponent 
in eq. (|7|) instead of i + 1 in eq. (||). Indeed summing these ‘double’ logarithmic singularities 
appears to lead to a singular behaviour / ~ x _A with A ~ 0.5 for qisrs and A q^s, and 
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A ^ 1 for the singlet densities (suggesting that the first moment of the singlet part of gi could 
be arbitrarily large). If one were to take these theoretical predictions seriously, the errors in 
the small x extrapolations considered below, particularly in the singlet channel, would have to 
be considered only as lower bounds. However the summation of ‘double’ logarithms is even less 
well founded theoretically than the summation of ‘single’ logarithms in the unpolarized singlet 
channel, and we believe that at present none of these results should be interpreted too literally 


Another important difference between the small x behaviour of unpolarized and polarized 
singlet distributions is that in the unpolarized case only the gluon anomalous dimension carries 
the leading singularity, and consequently the rise in the singlet quark distribution is driven 
directly by that of the gluon, while in the polarized case all the entries in the matrix of singlet 
anomalous dimensions are singular, and the polarized singlet quark and gluon distributions 
mix. It turns out that the leading eigenvector of small x evolution is then such that the singlet 
quark and gluon distributions have opposite sign, which means in practice that the singlet 
component of g\ is driven negative at small x and large Q 2 ||T6| . Contributions to first moments 
of g i from the small-x tail thus tend to become negative when Q 2 is sufficiently large. 


The purely perturbative asymptotic predictions eqs. (]6[.[7j) only hold when the input dis¬ 
tribution at the starting scale Q q is relatively nonsingular: if the singularity in the input is 
stronger than that generated perturbatively then the input will be essentially preserved by the 
perturbative evolution. The rise at small x will then be largely independent of Q 2 , rather than 
becoming steeper as Q 2 increases. If we take the starting scale in the crossover region between 
perturbative and nonperturbative dynamics, we can presumably take the small x behaviour of 
the input from Regge theory. For unpolarized distributions the input to the singlet distribu¬ 
tions (given by the pomeron trajectory) is then relatively flat, and indeed the dominance of 
the perturbative behaviour (]6j) is conErmed by F 2 data from HERA |TI| [2TJ] , while the input 
to the nonsinglet (given by the p — uj Reggeon trajectory) is singular, behaving as x -1 / 2 , so 
it is preserved by the evolution and is consistent with data from NMC and CCFR. For polar¬ 
ized distributions Regge theory suggests that the form of the input should be given by the A\ 
trajectory, and thus flat or even vanishing, behaving as x°-x 0,5 |I9] . 


In the following we will tentatively adopt a physical picture, inspired by the HERA results, 
in which we assume the validity of Regge behaviour at small x in the soft region (i.e. that 
at some input scale Qq ^ 1 GeV 2 the polarized densities are flat or vanishing) while at larger 
Q 2 the effect of NLO perturbative evolution is superimposed (giving a perturbative growth of 
the form eq. (0))- We will also allow steeper inputs in the nonsinglet sector, provided they 
are consistent with the partonic constraint \Aq N s( x , Q 2 ) \ ~ Qns { x , Q 2 ) : this limits the effective 
small x behaviour of A q^s to a singular growth ~ x -0 ' 5 . This picture turns out to be consistent 
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with the data, and gives a constraint on the allowed growth of \Aq N s(x, Q 2 ) \ in the unmeasured 
region which in turn limits the possible ambiguity on the Bjorken sum rule from the small x 
extrapolation. Unfortunately in the polarized singlet sector the corresponding partonic bound is 
much less restrictive and one has here to rely entirely on the validity of the NLO approximation 
and Regge behaviour at a low enough scale. 


2 Polarized Parton Densities from g 1 Data 


We now consider in detail the problem of extracting relevant physical quantities from the 
existing data. We devote particular attention to the study of the dependence of the results on 
the assumed functional form of the input parton distributions. For this purpose, we consider 
a variety of possible parameterizations of the input, we evolve these up to the values of x and 
Q 2 where data are available by solution of the evolution eqs. (Q) at NLO, and we determine 
the free parameters of the input by a best fit of gi(x, Q 2 ) eq. (||) to all the data of refs. j3|-[jrc[ 
with Q 2 > 1 GeV 2 . Experimental data for g\ are obtained from the experimentally measured 
asymmetries A\ using a parameterization of the measured unpolarized structure functions F 2 
f|33[| and R 


34)1 . P| consistently neglecting all higher twist corrections, and, for deuterium and 


helium targets, accounting for effects due to the nuclear wavefunction (but not Fermi motion 
or shadowing) by a simple multiplicative correction |35], 0]. Throughout this section we take 
a s {mz) = 0.118±0.005 and a 8 = 0.579±0.025 ]37| where a 8 is the SU(3) octet axial charge 


(in the proton, below charm threshold, Aq^sX) = Vns = § 9a + | as )- 

To begin with we parameterize the initial parton distributions at Q' 2 } — 1 GeV 2 according 
to the conventional form 


A/(x,Qjj) =N’ f rifX a f(l-x) f> t(l + 'y f x Sf ) (8) 

where Af denotes Aq^s, AT, or A g and A/} is a normalisation factor chosen such that the first 
moment of Af is equal to rjf. The signs of all parameters are left free, including the overall 
factors rjf (although the data always choose rjf to be positive). In table 1 we report the results 
from a fit of this class, denoted as fit A, which corresponds to: 

X = S g = 1, S NS = 0.75 (fixed), = 7 9 (fit A )- (9) 

It is particularly important to see to what extent the data fix the small x behaviour i.e. 
the exponents aj in eq. (H). Note that the rq data on protons and neutrons show a strong 
rise or fall at small x, while the deuteron data are much flatter. Thus the fitted nonsinglet 
quark densities at small x tend to show a sharp rise while the singlet quarks tend to remain 
fairly flat. Starting from a generic input set of densities of the class eq. (^) we can easily end 
up with A q^s considerably more singular than ~ x^ 0 5 , which we cannot accept since then 
the nonsinglet polarized distribution would eventually become larger than the unpolarized 

3 Although for the E154 data 0 the preliminary published values of g\ l are used directly since the experi¬ 
mental asymmetries are not yet available. 
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distribution. However a more careful analysis reveals that there is a strong correlation between 
aif and 5f. In the measured region one can easily push a^s closer to zero by decreasing 5 ns from 
unity without appreciable changes in the quality of the fit. For example for 5ns = 1, 0.75, 0.5 we 
find ajsrs ~ —0.8, —0.5, —0.3, respectively, with y 2 ~ 89, 89, 90. Thus we find that the existing 
data do not much constrain the behaviour of the nonsinglet at asymptotically small values 
of x: even within the simple functional form eq. (^|) one still has a considerable flexibility 
in the asymptotic behaviour as x —> 0. Still, it appears that, within this class of fits, the 
resulting nonsinglet quark density is rather more singular at Qp — 1 GeV 2 than we would 
expect from naive Regge theory, suggesting that Regge behaviour can only hold at significantly 
lower scales. The choice of parameters eq. (||) corresponds to a typical case which saturates the 
partouic constraint discussed at the end of the previous section. 


In order to discuss less singular inputs, we completely change the functional form of the 
input densities (while keeping the initial scale at Ql — 1 GeV 2 ). We thus choose an input in 
which the rise at small x is at most logarithmic (fit B): 


AS 
Acpvs 
A g 


A/s?/e^“ e (logl/x) /3s 

NnsVns (log 1/ x) aNS + (log 1/ x) i3ns (fit B), 

AfgTJg [(lOg 1/X)" 9 A 7 gX (lOg l/x)^ 


( 10 ) 


Note that logl/x (1 — x) near x = 1 , so that the 7 terms take care of the behaviour as 
x —> 1. The results we obtain from this fit are again reported in table 1. It is remarkable that 
we obtain a fit that is equally good in the measured region (as seen from the y 2 value) but 
with a quite different extrapolation in the unmeasured region. As a consequence, for example, 
the fitted value of ()a changes from gA = 1.17 ± 0.05 in fit A to qa — 1-23 ± 0.07 in fit B. The 
small x behavior of this fit is weaker than any power, and thus in particular compatible with 
the Regge prediction. 


Although the logs are reminiscent of QCD evolution the functional form of fit B might 
perhaps appear a little ad hoc. It is thus interesting to try to generate the logarithms in a more 
physical way, by perturbative evolution. In this spirit we consider another set of trials, where 
we start the QCD evolution at a very small scale, Q g = 0.3 GeV 2 , and fit a function of the 
form eq. (|8|) . The choice of such a low scale is simply used as a trick to generate an effective set 
of distributions at the value of Q 2 at which we begin to fit the data, i.e. Q 2 = 1 GeV 2 (data 
with lower Q 2 being still discarded), with the logs piled up in a way entirely consistent with 
perturbative evolution. In table 1 we report the results from a fit with 


7s = lg = 7 NS = 0 , Pg = 15 (fixed) (fit C). (11) 

In this class of fits, the large-x behaviour of the gluon distribution can hardly be determined 
by the fitting procedure; therefore, we fixed (3 g = 15 at Ql = 0.3 GeV 2 , because we checked 
that this choice approximately corresponds to a (1 — x) A behaviour of A g at Q 2 around 1 GeV 2 . 
Once more the quality of the fit in the measured region is unchanged. Comparing with the 
results of fit A, we see that by lowering the initial Qq scale all the exponents in the x a terms 
have become positive in qualitative agreement with the idea that naive Regge behaviour is 
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restored at a sufficiently low scale. Indeed a fit of comparable quality is obtained if we fix all 
exponents atf at Ql = 0.3 GeV 2 to the limiting value admitted by Regge theory, i.e. one half 
(fit D): 

a/= 0.5, 7s = 0, S g = 1, 6 ns = 1, (fit D). (12) 

The results of this fit are also shown in table 1. The y 2 is now slightly worse, but the physical 
results do not change much, especially in the nonsinglet sector (for example the central value of 
Qa is about the same in fits C and D). One could presumably optimize the choice of the initial 
scale Qq to make the agreement with Regge theory even better. 

In figs, la-c we display the best-fit g\ (fit B) for protons, neutrons and deuterons at the Q 2 
of the data. In figs. 2a-c we then compare the best-fit forms of g\ corresponding to fits A-D at 
Q 2 = 10 GeV 2 . The figures show that while the four fits are reasonably close together in the 
measured region (x ~ 0.003 — 0.03 up to x ~ 0.8), they become very different in the small x 
region. Note that at this Q 2 value g± has indeed become negative at small x in all cases. In 
fig. 3a-c we display the resulting polarized parton densities obtained from the fits at the same 
value of Q 2 . The behaviour at small x is quite different in each case. The fits C and D develop 
a particularly robust tail at small x in the singlet sector. It is the large positively polarized 
gluon that drives g\ negative at small x. 

In table 2 we report the values obtained by computing the first moments of g± by integra¬ 
tion of the four fits, both in the measured range of x and in the whole range at the ‘average’ 
values of Q 2 quoted by each experiment on protons, deuterons and neutrons. We see that 
while the truncated moments are remarkably close to each other the complete moments show a 
much wider spread. We also report the values of the truncated moments obtained by evolving 
the data to a common scale by means of the traditional (but unjustified) assumption that the 
asymmetries are scale independent and then summing over the bins, and those given by the 
experimental collaborations with their associated total errors. The latter two values should in 
principle coincide and only differ because of details in the way g\ is determined from the mea¬ 
sured asymmetries (such as the use of different parameterizations of the unpolarized structure 
function F 2 , or the inclusion of some higher twist corrections, as done by some experimental 
collaborations). The effect of the Q 2 dependence in the measured region is sizeable but smaller 
than the experimental error. Much larger is the indirect effect of scaling violations on the 
extrapolation at small x because of the larger scale dependence at small x. 


3 Phenomenological Implications 

We will now discuss the quantitative consequences that can be derived from the results of 
the previous section. In particular we will discuss the status of the Bjorken sum rule, the 
polarization of different quark flavors and of gluons in the proton, and the determination of a s . 



3.1 Testing the Bjorken Sum Rule 


One general way to test the Bjorken sum rule is to determine g A and the associated error from 
fitting the whole set of available data points. The error has many components to it: the error 
from the individual fit (based on a given functional form) due to the experimental errors on 
each data point, the ambiguity on the input functional form for the parton densities (which 
thus includes the error from the small x extrapolation), the error due to the uncertainty in 
the values of a s and as (we take a s (mz) = 0.118 ± 0.005 [flOH , and as = 0.579 ± 0.025 |37iP1), 


the error due to threshold uncertainties (estimated by varying the position of the thresholds 
around m q by a factor 2), and finally the error related to unknown higher perturbative orders 
(estimated by varying the renormalization and factorization scales /i 2 Rl gi 2 F by a factor of two 
in either direction). The error due to higher twist terms is potentially significant, because the 
data at low Q 2 (and thus at low x ) are essential for the computation of the first moments, but 
difficult to estimate reliably: here we estimate it (see sect. 3.4 below) from the shift in the value 
of gA extracted from the Bjorken sum rule due to the inclusion of a higher twist correction to 


it according to QCD sum rule and renormalon estimates |39 


Based on the fits shown in table 1, and similar ones used to estimate theoretical errors, we 


find 


g A = 1.19 ± 0.05(exp) ± 0.07(th) = 1.19 ± 0.09, 


(13) 


where the central value is obtained as the average between the maximum and minimum values 
of table 1, and the various contributions to the theoretical error are listed in table 3. The fitted 


value is to be compared with the direct measurement g A = 1.257 ± 0.003 |40| from /1-decay. 
Thus we find that the Bjorken sum rule is confirmed to within one standard deviation but still 
with an accuracy of only about 8%. For the chosen value of a s the fitted value of g A is below 
the experimental value; it would be increased by using a larger value of a s , or perhaps by the 
inclusion of additional corrections beyond NLO. 


3.2 Singlet First Moments 


Similarly in the singlet sector one obtains from the data values for r) q = A£(l) (the conserved 
polarized singlet quark density), for ?j g = Ag(l,lGeV 2 ) (the first moment of the polarized 
gluon density evaluated at Q 2 = 1 GeV 2 ), and for a 0 (10 GeV 2 ) (the non conserved singlet axial 
charge defined implicitly from the singlet part of gi by eq. (§)). This latter quantity approaches 
a finite limit at infinite Q 2 because the corresponding anomalous dimension starts at two loops, 
and within the present accuracy ao(10 GeV 2 ) is equivalent to ao(oo). The values for these three 
quantities as obtained from our representative fits are reported in table 1. We then studied in 
detail, following ref. [17]], the theoretical errors from the various different sources: results for 
these are listed in table 3. The error due to higher twists is not included since it cannot be 


4 The sensitivity to as is in practice very weak: we have checked that variations of up to 30% have little 
effect on our results. Such deviations from the value of ref. Q could be caused by SU(3) violations in hyperon 
semileptonic decays Q. 
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reliably estimated but, on the basis of the corresponding error on ()a, it is expected to scarcely 
affect the total errors which in the singlet sector are rather large (for reasons discussed in detail 
in refs. [[R], [L7]]). In comparison to ref. 13 we have a more reliable determination of the error 
related to the fitting procedure, and a somewhat stronger sensitivity to higher order corrections. 
We then find 


A£(l) 

As(l,lGeV 2 ) 


a 0 (oo) 


0.45 ± 0.04 (exp) ± 0.08 (th) = 0.45 ± 0.09, 
1.6 ± 0.4 (exp) ± 0.8 (th) = 1.6 ± 0.9, 

0.10 ±0.05 (exp)lg;^ (th) = 0.10tg;ft, 


(14) 


in excellent agreement with the corresponding results of ref. (T^J. The physical implications 
are thus the same. The parameter ao(oo) measures the degree of ‘spin crisis’: the singlet axial 
charge of the nucleon is still compatible with zero as it was at the beginning of the story 
Note that with the naive Regge extrapolation at small x the experimental result for the 
axial singlet charge would be significantly larger, with a much smaller error: for example in 
ref. 0 a value ao(oo) = 0.33 ± 0.04 was quoted. There is also some evidence (around two 


standard deviations) for a positive gluon polarization in the nucleon (increasing with Q 2 as 
1 /a s (Q 2 )). The amount of gluon polarization is large enough to allow the first moment AE(1) 
of the conserved singlet quark density to be within one standard deviation of as ~ 0.58, which 
in the absence of all SU(3) and chiral symmetry breaking effects, could be identified with 
the constituent spin fraction J|!J. This can be seen as a direct confirmation of the physical 
explanation of the ‘spin crisis’, advocated in refs. as due to the anomaly and well 

described in terms of the QCD parton model (for more general possibilities, see refs. ||43|| ). 


3.3 Determination of a s 

The above discussion on (Ja makes it clear that the determination of a s from the Bjorken 
integral is adversely affected by the increased ambiguity from the small x extrapolation that 
follows from the demise of the naive Regge behaviour at small x. Fortunately we find that a s 
can be determined directly from the available data without extrapolation in the small x region 
if the totality of the data is taken into account and not just the Bjorken integral. The value of 
a s is then determined by the strong scaling violations needed to accommodate the difference 
between the data at small Q 2 from the SLAC experiments and those at larger Q 2 from the 
SMC in the common range of x. While Ag(l, 1 GeV 2 ) is mainly fixed by the proton data, a s 
is determined by the difference between proton and neutron. To show this we repeated the 
fits A-D but fixing qa to its experimental value and instead fitting a s (in this case, fit A was 
performed with j3 g fixed at the value /3 g — 4). In all cases the central value was found to be close 
to a s {mz) = 0.120. Since the different fits differ considerably in the unmeasured region, this 
shows that it is the behaviour in the measured region that matters. In addition, the Bjorken 
integral is appreciably different in the different fits A-D when a s is kept fixed (and the value 
of Ag(l,lGeV 2 ) even more so) but this difference does not affect the fitted a s very much. 
Furthermore, the resolution of the discrepancy between the fitted value of qa eq. 0) and its 
experimental value would require a much larger increase of a s if only the Bjorken integral were 
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relevant for fixing a s . These results show that a s is much better constrained by the overall 
pattern of scaling violations than by the Bjorken integral alone. 


The theoretical uncertainties that affect this determination of a s are listed in table 3. The 
main source of uncertainty originates from higher order and higher twist contributions. Higher 
order corrections are estimated by varying the renormalization and factorization scales /i 2 Rl [i 2 F 
around Q 2 . The values of table 3 are obtained varying /j R and /jl 2 f between Q 2 / 2 and 4 Q 2 . A 
wider scale variation at the lower edge is made impossible, without leaving the perturbative 
domain, by the presence of data with values of Q 2 close to 1 GeV 2 . The higher twist error is 
again determined from the corresponding shift in the value of a s extracted from the Bjorken 
sum rule (sect. 3.4 below). We can cross-check this estimate of higher twist terms by repeating 
the fit while excluding all data with Q q < 2 GeV 2 . We find that the value of a s does not 
change, thus suggesting that higher twist effects are indeed quite small. 

We thus obtain finally 

a,(m z ) = 0.12<V”;°“(expU° : ““(th) = 0.120 ±g;8Jg. (15) 


This reasonably good determination of a s {mz) could still be improved with better data: it is 
important to notice that without the very recent neutron data m the experimental error would 
be twice as large. We see no reason why it could not be as good as the determination from 
unpolarized data if more data were added and the experimental errors consequently reduced. 
However the theoretical error is already the dominant one; it could also be reduced by more 
data because it is made particularly large by the small Q 2 values of the neutron measurements. 


3-4 Determining a s from the Bjorken sum rule 


Even though the full set of anomalous dimensions and coefficient functions are known at NLO, 
for the nonsinglet first moment the anomalous dimension vanishes and the coefficient function 
is known to order cc^, i.e. at NNNLO J44|. One would like to take advantage of this situation 
by directly determining a s from the Bjorken sum rule p3|. This however requires first a 


determination of the Bjorken integral from the data. Since the latter can only be accomplished 
at NLO the procedure is only advantageous to the extent that this extraction can be done in 
a way that minimizes the necessary theoretical input. 

This can be done following a procedure close to that used in the experimental papers: first 
one computes directly from the data the contribution to the first moment from the experimen¬ 
tally accessible x range (applying a NLO correction to account for evolution of the data at a 
common scale for each experiment), then one adds to it an extrapolation over the unobserved 
x-range. The resulting values of the Erst moments are finally analysed using the NNNLO result 
for the Bjorken sum rule in order to extract a s . The procedure can be justified if the Q 2 span of 
the data around the average value in each experiment is small in comparison to the difference in 
average Q 2 between experiments. This is indeed the case for most of the data i-G3 , with the 
exception of the SMC data at very large or very small x, which however give a relatively small 
contribution to the first moments. Although it is generally inconsistent to correct the data 
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by using scaling violations at NLO and then to pretend one has a NNNLO accuracy for first 
moments, in practice this might not be a cause for real concern since the corrections to the data 
on the truncated first moments from the NLO scaling violations are not large in comparison 
with the experimental errors and those from the extrapolation in the unmeasured region (see 
table 2). 

We thus first determine for each experiment the value of the first moment of gi in the 
measured region, evaluated at the average scale ( Q 2 ) quoted for that particular experiment, by 
simply summing the experimental determinations gf over the bins. Corrected values of g\ in 
order to account for evolution at a common scale are determined according to || 

9i(x, (Q 2 )) ~ g{{x, Q 2 ) + A(x, Q 2 , (Q 2 )), (16) 

A(x,Q 2 ,(Q 2 )) = g[(x, (Q 2 )) - g((x,Q 2 ), (17) 

where g{(x,Q 2 ) is a best-fit, which we take from our sets of fits A-D. The fits with g^ left 

as a free parameter must be used here, otherwise the determination of the nonsinglet first 
moment would be circular. Because all these fits give very similar results in the measured 

region (compare table 2), we can assume the error in the evolution procedure to be negligible 

and take the error in this contribution to the moment to be that given by the experimental 
collaborations. We then complete the moments in the unmeasured region by again using the 
set of fits A-D, which differ widely in the small x region, but are all consistent both with 
the experimental data and well understood theoretical principles. Specifically, we take (for 
each experiment) the average of the difference between the full and truncated integrals for the 
four fits as an estimate of the extrapolation. The associated uncertainty is then taken to be 
given by their spread as one varies the functional form of the fit, as well as renormalization 
and factorization scales, and the value of a s (variations of the thresholds and the value of a§ 
have no significant effect). The measured and extrapolated first moments are then added and 
the respective errors added in quadrature (see table 4). This is then the best estimate of the 
experimental determination of the first moments and their associated errors. 

In order to determine the isotriplet combinations relevant for the Bjorken sum rule we have 
repeated the same analysis directly at the level of the isotriplet combinations of first moments, 
which we can determine separately at 10 GeV 2 for the SMC experiment, and at 3 GeV 2 for all 
of the SLAC experiments (by averaging the neutron data). The error due to the extrapolation 
is then only sensitive to the uncertainty in the nonsinglet contribution to g\ at small x and 
thus much smaller than the sum in quadrature of the errors in individual first moments: in 
practice only the variation of the functional form of the fit contributes, because the nonsinglet 
extrapolation is essentially insensitive to the value of a s and the choice of renormalization and 
factorization scales. The result, listed in table 4, shows that whereas for the SMC data the 
dominant uncertainty is in the measured region, for the SLAC data the uncertainty in the 
extrapolation is comparable to it. 

From the isotriplet first moments of table 4 we can then determine a s (mz ) by fixing to 
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SMC only, 
SLAC only, 
all data, 


(18) 


its value qa = 1.2573 ± 0.0028, eq. ([13]). We find 

f O.llUHii (exp) 
«.N = < 0-118^8:02° ( ex p) 
( 0.118±8:3m (exp) 


where for the SMC only case the lower error should be understood as meaning that the un¬ 
certainty is 100% or more. The error in eq. (|H|) is the purely experimental one due to the 
error on the isotriplet first moment of table 4 . This error is now so large that it dominates 
any theoretical errors. For instance, we can estimate a higher twist error by adding to the 
Bjorken integral a twist four term, with a coefficient estimated from sum rule and renormalon 
calculations []D| . The value of a s of eq. (|T%| ) is then reduced by about 0.004. We used this as an 
estimate of the error from higher twist terms in the determination of a s eq. © (see table 3). 


As a cross-check, we can instead fix a s to its known value a s (mz) = 0.118 ± 0.005 and 
extract qa from the sum rule: we then get g\ = 1-27 ± 0.11 ± 0.05 ± 0.03 where the first error 
is experimental, the second is due to the error on a s , and the third is the higher twist error, 
estimated as above (and used in table 3). These errors are comparable to those found by the 
direct determination of qa eq. (p~3|) , but the experimental error is somewhat larger because only 
a subset of the information contained in the data is now being used. 


We conclude that the result for a s obtained from this procedure is compatible with that from 
the direct determination eq. (0). but that the error is now considerably larger, in agreement 
with our previous observation that a s is better constrained by the data in the measured range 
than by the Bjorken integral alone. The value of ot s (mz ) is very close to that found in ref. (23[ 


but the error is larger because of the error from the small x extrapolation which we no longer 
assume to follow naive Regge expectations at ( Q 2 ). A direct determination of a s from the 
Bjorken sum rule could however be competitive if data with a wider kinematic coverage in x 
were available. 


4 Conclusions 


We have performed a global analysis of all data on the polarized structure function g\ for proton, 
deuteron and neutron targets, including data which have only recently become available, using 
NLO perturbative QCD. We took care to properly estimate the effects of perturbative evolution 
and the uncertainties in the small x extrapolation when computing first moments. We showed 
that the data confirm the Bjorken sum rule to within one standard deviation at the 8% level, 
and indicate a gluon polarization in the nucleon which is non-zero at the level of two standard 
deviations. We showed further that the data now provide a reasonably accurate determination 
of a s , consistent with the global average. However we also showed that the Bjorken integral is 
unfortunately not yet sufficiently well determined to admit a competitive determination of a s 
using the Bjorken sum rule. 

Our determination of a s could be improved by more accurate data on polarization asym¬ 
metries in the fixed target region, or by an independent (semi-inclusive) determination of the 
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polarized gluon distribution. However the theoretical error from higher order corrections is 
already the dominant one. The error in our determination of first moments is instead domi¬ 
nated by uncertainties in the small x region, and we suspect that these uncertainties will only 
be significantly reduced by the measurement of g\ at lower x and higher Q 2 , which would be 


possible at a polarized colliding beam experiment at HERA [(45 


Acknowledgement: We thank P. Bosted and A. Deshpande for useful information on the 
experimental data. 
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Parameters 

A 

B 

C 

D 

d.o.f. 

114-11 

114-11 

114-8 

114-8 

Ql/GeV 2 

1 

1 

0.3 

0.3 

V* 

0.408 ± 0.041 

0.410 ±0.039 

0.422 ±0.026 

0.492 ± 0.036 

as 

0.741 ±0.353 

1.710 ±0.416 

2.600 ±0.964 

0.5 (fixed) 

Pt, 

3.105 ± 1.049 

2.735 ± 0.466 

3.359 ± 1.210 

1.039 ±0.241 

7e 

0.185 ±2.496 

0 (fixed) 

0 (fixed) 

0 (fixed) 

Vg 

1.068 ±0.403 

1.032 ±0.330 

0.479 ± 0.095 

0.650 ±0.104 

a g 

-0.597 ±0.286 

2.970 ±0.611 

0.217 ±0.319 

0.5 (fixed) 

Pg 

0.831 ±2.322 

1.286 ±0.895 

15 (fixed) 

13.19 ±9.57 

7 9 

0.185 ±2.496 

21.2 ±22.5 

0 (fixed) 

-0.548 ±9.139 

9a 

1.168 ±0.052 

1.234 ±0.066 

1.146 ±0.038 

1.141 ±0.036 

as 

0.579 (fixed) 

0.579 (fixed) 

0.579 (fixed) 

0.579 (fixed) 

OiNS 

-0.537 ±0.057 

1.656 ±0.166 

0.765 ±0.228 

0.5 (fixed) 

Pns 

2.503 ± 0.274 

5.320 ±0.251 

2.087 ±0.448 

2.622 ±0.410 

Ins 

17.46 ± 8.43 

-0.229 ±0.105 

0 (fixed) 

5.024 ±4.647 

x 2 

83.7 

83.0 

83.8 

90.9 

X 2 / d.o.f. 

0.813 

0.806 

0.790 

0.858 

As(l,lGeV 2 ) 

1.07 ±0.40 

1.03 ±0.33 

1.61 ±0.32 

2.08 ±0.34 

a 0 (10GeV 2 ) 

0.15 ±0.07 

0.16 ±0.05 

0.05 ±0.04 

0.02 ±0.04 


Table 1: Results of fits A-D described in the text 



SMC : p 

E143 : p 

SMC : d 

E143 : d 

E142 : n 

E154 : n 

(Q 2 )/GeV z 

10 

3 

10 

3 

2 

5 

Meas. Range: Exp. 

0.1310 

0.1200 

0.0379 

0.0400 

-0.0280 

-0.0370 

Exp. Error 

±0.0156 

±0.0089 

±0.0079 

±0.0050 

±0.0085 

±0.0108 

Ai ind. Q 2 

0.1350 

0.1059 

0.0442 

0.0398 

-0.0293 

-0.0364 

Meas. Range: A 

0.1256 

0.1070 

0.0430 

0.0378 

-0.0328 

-0.0338 

Meas. Range: B 

0.1268 

0.1061 

0.0425 

0.0373 

-0.0326 

-0.0342 

Meas. Range: C 

0.1280 

0.1084 

0.0462 

0.0401 

-0.0289 

-0.0316 

Meas. Range: D 

0.1308 

0.1083 

0.0492 

0.0402 

-0.0306 

-0.0308 

Full Range: A 

0.1171 

0.1137 

0.0275 

0.0265 

-0.0600 

-0.0614 

Full Range: B 

0.1232 

0.1197 

0.0286 

0.0276 

-0.0637 

-0.0654 

Full Range: C 

0.1039 

0.0999 

0.0160 

0.0144 

-0.0707 

-0.0716 

Full Range: D 

0.0990 

0.0946 

0.0115 

0.0094 

-0.0754 

-0.0759 


Table 2: Determination of the first moment Ti((Q 2 )) eq. (||). For each experiment we display 
the average value of Q 2 and the contribution to the first moments from the measured range of x, 
as given, first, by the experimental collaborations, with the corresponding total (statistical and 
systematic) error, then by summing over experimental bins while evolving the data assuming 
scale independent asymmetries, and finally as obtained from integration of the fits A-D. In the 
last four rows the complete first moments obtained from the fits A-D are shown. 
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9A 

AS 

A g 

a o 

OL s 

experimental 

±0.05 

±0.04 

±0.4 

±0.05 

+0.004 

-0.005 

fitting 

±0.05 

±0.05 

±0.5 

±0.07 

±0.001 

a s & as 

±0.03 

±0.01 

±0.2 

±0.02 

±0.000 

thresholds 

±0.02 

±0.05 

±0.1 

±0.01 

±0.003 

higher orders 

±0.03 

±0.04 

±0.6 

+0.15 

-0.07 

+0.007 

-0.004 

higher twists 

±0.03 

- 

- 

- 

±0.004 

theoretical 

±0.07 

±0.08 

±0.8 

+0.1V 

-0.010 

+0.009 

-0.006 


Table 3: Contributions to the errors in the determination of the quantities Qa, AE(1), 
Ag(l, lGeV 2 ), ao(oo) and a s (mz) from the fits described in the text. 



m 

Meas. Range 

Full Range 

SMC : p 

10 

0.133 ±0.016 

0.116 ±0.022 

SMC : d 

10 

0.045 ± 0.008 

0.021 ±0.016 

SMC : 1=1 

10 

0.176 ±0.036 

0.191 ±0.037 

E143 : p 

3 

0.108 ±0.009 

0.107 ±0.017 

E143 : d 

3 

0.040 ± 0.005 

0.021 ±0.014 

E142 : n 

2 

-0.032 ±0.009 

-0.068 ±0.015 

E154 : n 

5 

-0.034 ±0.011 

-0.070 ±0.015 

SLAC : 1=1 

3 

0.139 ±0.013 

0.177 ±0.018 


Table 4: Best estimates and errors for the truncated and the complete first moments for the 
target and the average Q 2 that correspond to each experiment. The isotriplet first moments 
for the SMC experiment, and for all the SLAC experiments combined are also given. 
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0.00 
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dashed: q 2 = q 2 (x) (SMC) 
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dotdashed: q 2 =l GeV 2 
dotted: q 2 =10 GeV 2 
x : SMC 
O: E143 







Figure 1: Plots of xg±(x, Q 2 ) for fit B for(a) proton (b) deuterium (c) and neutron 
targets. The curves correspond to Q 2 = 1 GeV 2 (dotdashed), Q 2 = 10 GeV 2 (dotted), 
and Q 2 {x) of the various experiments: Elf3 (solid) and SMC (dashed) for figs, a-b 
and Elf 2 (solid) and Elf 5 (dashed) for fig. c. The data points with total errors are 
also shown. 
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